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PACS 03 . 75 . Be - Atom and neutron optics 
PACS 03 . 75 . Pp - Atom lasers 

PACS 67 . 85 . Jk - Other Bose-Einstein condensation phenomena 

Abstract. - Given a quantum state at some instant of time t, the underlying system Hamiltonian 
can not only predict how the state will evolve, but also the history of the state prior to t. Thereby, 
in order to have a directed motion, like in a diode, some sort of irreversibility must be considered. 
For the atom diode, this has been achieved by spontaneous decay of excited atomic levels. For an 
atomic condensate, it is clear, however, that such decay will cause both heating and decoherence of 
the condensate. To overcome this comphcation we introduce a different setup where the dissipation 
does not act directly on the atoms. The excited atoms are stimulatedly driven back to the ground 
state by exciting a cavity mode, which in return decays to the vacuum via photon losses. The 
efficiency of the method utilizing experimental parameters is shown to be almost perfect within 
large parameter regimes. 



Introduction. — With the advances in atomic cool- 
ing [2] and in the construction of atomic chips [3], it is 
nowadays possible to coherently guide and manipulate 
atomic matter waves [31 [5]. Retaining the inherent co- 
herence of the matter waves is a necessary condition for 
successful applications. This becomes maybe especially 
important in schemes of quantum information processing 
with neutral atoms 6 and for quantum memories 7 . The 
building blocks in these atom optics models are formed by 
light-matter interaction in which the light fields constitute 
effective matter wave potentials. By carefully designing 
the light-matter interaction one can achieve analogs of for 
example lenses [Sj , transistors and amplifiers [9j , and atom 
diodes [10Hl2]. 

In this paper we consider an atom diode - a device that 
allows atoms in their ground state to pass in one direction 
but not in the opposite direction - where decoherence due 
to irreversible losses are greatly suppressed. One idea of 
using atomic Bose-Einstein condensates (BEC) in atom 
optics setups is the fact that effective time-scales (often 
scaling inversely with the number of atoms) can be con- 
siderably decreased, and moreover, the many atomic inter- 
nal degrees of freedom of the BEC permitting storing of a 
large amount of quantum information. However, conden- 
sation relies on coherence and it is thereby desirable to 
minimize the amount of decoherence in the atom diode. 
Our idea is based on the one of ref. [10] (see also [13]), 
but we remove atomic spontaneous emission and replace 



it by Raman coupling the atoms to a high-Q optical cavity. 
Spontaneous emission, as utilized in ref. [10 , will induce 
both heating of the condensate as well as loss of inherent 
coherence, while in our scheme the atoms follows a dark 
state during the Raman process and the atomic motion 
is, within the adiabatic regime, constant. In addition, the 
irreversible process occurs between the cavity field and 
its surrounding photon reservoir, and consequently there 
are no direct decoherence effects on the condensate. Em- 
ploying experimentally relevant parameters, we find that 
the diode efficiency is close to perfect within a mean-field 
treatment. Non-adiabatic contributions, arising either due 
to high atomic velocities or strong non-linearities, result in 
atom- field entanglement. These correlations decay thanks 
to the reservoir-induced collapse of the cavity field, which 
therefore improves the overall efficiency. However, tunnel- 
ing resonances may cause non-zero transmission of back 
scattering atoms. This, on the other hand, can be con- 
trolled, and consequently suppressed, by adjusting the ex- 
ternal lasers. 

Model system. — We consider a BEC of ®^Rb-atoms 
propagating along the positive direction in a waveguide 
aligned in the x-direction, i.e. the dynamics in the trans- 
verse directions can be frozen out. However, the trans- 
verse confinement is still not large enough to cause neither 
confinement-induced resonances [14] nor collapse of the 
condensate [15]. Two atomic ground state Zeeman levels 
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|1) and |3) (typically taken from the 55*1/2 ^ = 1 Zeeman 
manifold) are Raman-coupled to an excited state |2). Two 
spatially separated classical lasers drive the 2 o 3 transi- 
tion and another classical laser the 1 •(-> 2 transition. In 
addition, the atomic states |1) and |2) are also coupled via 
a mode of an optical Fabry-Perot resonator. The atomic- 
laser couplings are schematically presented the level dia- 
gram of fig. [1] (a). All coupling fields propagate in the 
transverse direction and thereby possess Gaussian profiles 
in the longitudinal x-direction. Furthermore, they are all 
detuned from their respective atomic transitions by A. Fi- 
nally, the state |1) is Stark-shifted by yet another classical 
laser. The positions of lasers along the waveguide and 
their relative strengths are displayed in fig. [1] (b). Thus, 
the condensate traverses the pulses in the order; a Stokes 
(s) pulse coupling states 2 and 3, a pump (p) pulse cou- 
pling 1 and 2, the Stark pulse acting on state 1, followed 
by another Stokes pulse formed by the cavity and coupling 
states 1 and 2, and finally a pump pulse between 2 and 3. 
The amplitude of the Stark-shift pulse W has been taken 
twice as big compared to those of the Raman Stokes and 
pump pulses i}s,p and Gs,p, 
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Fig. 1: Atomic level configuration (a) and positions of the 
lasers along the waveguide (b). The two pulses fip and Gs 
couple the ground state |1) to the excited state |2), while the 
remaining two pulses fip and Gs couple instead the atomic 
states 1 2) and |3). The detuning A prevents population of the 
excited state for atoms propagating in the negative direction. 
The Stark-shift pulse W, acting solely on the atomic |l)-state, 
has twice the amplitude of the Raman pulses, while all the 
widths are the same. 
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G fj{x) ^ fiVtoe , p^p^s. 

The pulse positions Xg , Xp and ys , j/p are chosen such that 
population is adiabatically swapped between the two Zee- 
man ground state levels |1) and |3) via the stimulated Ra- 
man adiabatic passage (STIRAP) process [TB]. To achieve 
a robust population transfer, the two Raman pulses should 
occur in a counterintuitive way and overlap |17| . The ad- 
vantages of an adiabatic scheme, provided the evolution is 
adiabatic, is that it is insensitive to parameter strengths 
as well as effective interaction times. Yet another advan- 
tage using a STIRAP scheme is that the evolution follows 
a dark state tpdix) such that the condensate velocity is 
ideally not afi^ected by the couplings. As long as we con- 
sider ultracold atoms, the center Stark-shift pulse hinders 
atoms in state [1) to pass. However, atoms entering from 
the left in state |1) will first be STIRAP transfered into 
|3)-state atoms before arriving at the center pulse, and 
hence pass the barrier unaffected. The second STIRAP 
sequence reverses the population swapping, but this time 
the Stokes pulse consists of a high-Q optical cavity mode 
such that each atom passing it leaves behind a photon. On 
the time-scale of the atomic motion, the photons rapidly 
escape the cavity leaving it in vacuum. Therefore, |1) 
atoms coming from the right cannot be transfered into the 
|3) states since the cavity is empty from photons. Atoms 
arriving from the right in the |3) state, on the other hand, 
couple to the excited |2) atomic state via the Pump pulse 
Gp{x). However, this transition is highly detuned due to 



the non-zero A and hence the population transfer greatly 
suppressed. The atoms in this state, nevertheless, experi- 
ence an effective repulsive Stark shift potential upon which 
they are reflected. Consequently, |1) atoms entering from 
the left are reflected by the W{x) Stark shift pulse, while 
3) atoms are reflected by the detuned Gp{x) pulse. The 
amount of thermally excited atoms in the |2} state is as- 
sumed vanishingly small. 



For the F ~ 1 manifold of ^^Rb atoms, the s-wave spin- 
dependent interactions is a factor 200 smaller than the 
spin-independent interactions [18]. In the present situa- 
tion, the effect of interactions plays a minor role and we 
can safely neglect the spin-dependent interactions. Denot- 
ing the s-wave scattering length as as, the atomic mass as 
TO, and the transverse condensate widths At, the effec- 
tive atom-atom interaction reads U — 2h^as/mA^ [15] . 
Within a mean-field approach, the Gross-Pitaevskii equa- 
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tion in the rotating wave approximation, 
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render the time-evolution of the order parameter. Here, 
the spin operators cTij — |i)0|, N is the total number 
of atoms, a and represent the photon annihilation and 
creation operators of the cavity mode, and the spinor order 
parameter 

ipw{x,t) 
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where \ipij{x,t)\'^ characterizes the atomic density of level 
\i) and with an empty j = 0, or populated cavity j = I. 
We have normalized the order parameter to unity, and 
thereby the appearance of the factor N in the interaction 
term. Note that the coupling schemes of the Stokes and 
the pump pulses have been swapped between the two Ra- 
man processes in order to reverse the population transfer 
among the Zeeman levels |1) and |3). 

We have treated the condensate at a mean-field level, 
while the cavity field is at a quantum level. We could re- 
place the operators a and by the coherent amplitudes 
a and a* respectively in the equation above. However, 
when treating cavity losses within a master equation ap- 
proach we need to keep the quantum nature of the cavity 
field 0. We take a and as the regular boson operators, 
but it is understood that while the condensate traverse 
the cavity field its true amplitude is proportional to the 
atom number N. In our model, the cavity field ampli- 
tude is independent of the atomic number A'' and instead 
jap = 1. The effective model that we consider for the 
cavity field should not be a sever shorting since the cavity 
field decay is independent on the field amplitude, namely 
a{t) = a(0)e~''* where k is the decay rate. Moreover, The 
crucial aspect is not in the actual cavity decay as long as 
it is fast compared to the time-scale set by the condensate 
motion. 



^Solving the effective master equation takes into account for 
reservoir-induced quantum noise, which typically is absent in a full 
mean-field approach. 



As pointed out in the introduction, for a diode we should 
have a directed motion of the condensate, i.e. if we reverse 
time we should not be able to propagate to negative x's 
provided the condensate is located at large positive x's. 
Despite rendering unitary evolution, the Gross-Pitaevskii 
equation is inherently non-linear and time-reversal may 
cause hysteresis effects. This appears only for very high 
atomic densities or strong atom-atom interactions, which 
does not apply to our model. The irreversible dynamics 
arises due to cavity losses, which is here modeled by a 
master equation. For the density operator initially in the 
state p = 'i'^x, 0)'i^{x, 0), and for coupling to a zero tem- 
perature photon bath, the equations of motion takes the 
form [19] 
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Losses of the atomic Zeeman levels have been neglected 
as the excited |2)-state is negligible populated within the 
adiabatic regime that we consider. 

Results. — The numerics will be for a ^^Rb conden- 
sate as for the BEC-cavity experiments of Ref. [20] . More 
precisely, we take a condensate of iV = 100 000 atoms, the 
transverse widths At = 3 pm, scattering length as = 5.77 
nm, atom-field couplings g ^ Hq = 2tt x 10.9 MHz, detun- 
ing A = 27r X 55 MHz, cavity mode waist w = 15 pm, and 
cavity loss rate k = 2tt x 1.3 MHz. Note that via eq. 
the atom-field coupling and the mode waist set the ampli- 
tudes and waists for the classical lasers as well. The pulse 
locations are Xp = —ya — —130 /im and Xg — —Up — —160 
pTn. As an initial state for a condensate propagating in the 
positive direction we take a Gaussian located at xq = —260 
/im, with a longitudinal width A; — 10 pva, an initial ve- 
locity in the positive a;-direction of uo = 5 cm/s [21j . and 
only the |l)-component populated 
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The actual shape of the condensate has little infiuence on 
the protocol, i.e. we could equally well have chosen some 
sort of Thomas-Fermi profile. The cavity mode is initially 
in vacuum. For verifying the diode efficiency we also con- 
sider condensates propagating in the negative direction. 
Then the initial state is the same but with xq —xq and 
vq — )• —vq. we also check for negatively propagating atoms 
initially in the |3) state and with the same motional wave 
packet. 

We use a Monte-Carlo method combined with 
the split-operator method [23] to solve the master equa- 
tion (jll). The initial state is propagated with the non- 
Hermitian Hamiltonian H = Hop — ^a^a. This propa- 
gation is performed on some time- grid with steps dt. For 
this we utilize the split-operator method, and chose dt 
small enough such that convergence is guaranteed. At 
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Fig. 2: The plots (a), (c), and (e) show the populations Pi{t) of 
the different Zeeman-levels |1) (solid curve), |2) (dotted curve), 
and 1 3) (dashed curve). The condensate velocity v{t) is de- 
picted in (b), (d), and (f). In the first two plots we give the 
result for atoms propagating from left to right. In this case, 
the first STIRAP-process results in almost perfect population 
transfer, indicating the adiabatic following of the dark state. 
The velocity fluctuations in the second STIRAP-process is due 
to transitions out of the dark state induced by the cavity de- 
cay. The results of atoms initially in state |1) propagating in 
the negative direction (from right to left) are shown in (c) and 
(d). Here it is clear how the atoms remain in the same Zeeman- 
level and how they are reflected by the W{x) pulse. Finally, 
in (e) and (f) we present the same but for atoms initially in 
the 1 3) state. This time the atoms are reflected already at the 
Gp{x) pulse. The parameters are given in the text. 

each time-step we randomly simulate a quantum-jump de- 
scribing loss of a photon. One such simulation render a 
single quantum-trajectory. Repeating this procedure and 
average over all trajectories results in the desired solution. 
We focus on the population of Zeeman levels 



the velocity 



dxit) 



where x(t) = Tr xp{t) , and the reduced densities 



p,{t) = Trj^, p{t) 



(6) 
(7) 

(8) 



In this last equation, trace is over the cavity field and the 
remaining two Zeeman levels' degrees of freedom, i.e. Pi{t) 



represent the atomic density of Zeeman level Since 
the atomic velocity should ideally not couple to any other 
degrees-of-freedom, the velocity is an indirect measure of 
adiabaticity. More precisely, the initial state populates a 
dark state with the corresponding adiabatic potential [53] 
being strictly zero, and as long as the dynamics is adi- 
abatic the velocity should stay constant. On the other 
hand, if any of the bright states become populated due 
to non-adiabaticity it will render a velocity change since 
their corresponding adiabatic potentials are non-zero and 
thereby generates an effective force. The bright states nor- 
mally contain as well population of the Zeeman state |1) 
which will be reflected upon the center pulse W{x). In 
addition to the above quantities, we will also calculate the 
total transmission probability T in order to have a second 
verification of the directed motion of the condensate. 

Figured] presents various numerical results of the popu- 
lations and the velocities, obtained using the above men- 
tioned parameters and averaging over 500 Monte Carlo 
trajectories H. The first two plots, (a) and (b), give the 
results for atoms propagating from left to right. Within 2 
% accuracy, the final population in the target state [1) and 
the transmission are perfect. Except for the second Ra- 
man process, the velocity stays roughly constant signaling 
adiabatic evolution. During the second Raman process, 
as photons decay, the system cannot follow the dark state 
ipd{x) adiabatically. Consequently, since the decay is not 
an adiabatic process, the velocity changes within this re- 
gion. However, once leaving the cavity interaction zone, 
the initial velocity vq is restored, as will be discussed fur- 
ther at the end of this section. 

In the following two plots, (c) and (d), we display the 
same but for a condensate in the internal |l)-state prop- 
agating in the negative direction. This time, since the 
cavity is empty the condensate traverses the first cou- 
pling pulses unaffected and reflect upon the Stark-shift 
pulse W{x). Finally in (e) and (f), the results for neg- 
atively propagating condensates initially in the |3)-state 
are shown. Due to the large detuning, A w h5flo, approx- 
imately no population transfer between the atomic states 
2) and |3} takes place. Instead, the Gp{x) pulse acts as 
a Stark-sfift and the condensate gets reflected by it. For 
the examples of (c)-(f), the transmitivity T is of the or- 
der 10~^. If the atomic velocity is increased, atoms will 
at some point traverse the Stark-shift pulses. For atoms 
propagating from right to left, we have numerically seen 
that for velocities vq < —50 m/s the atoms in the |1) state 
begin to pass the reflective W{x) potential. At this ve- 
locity, the transmission is about 2 %. For the employed 
system parameters, we do not find any tunneling reso- 
nances [25 within the velocity interval 5-50 cm/s. Beyond 
velocities vq < 50 cm/s, the |l)-state atoms have enough 
kinetic energy to pass the central barrier. The reflective 
pulse Gp{x) for |3) atoms is not as strongly reflective as 



^For backward propagation, the cavity stays in vacuum and no 
quantum jumps occur. Consequently a single trajectory captures 
the evolution in these cases. 
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Fig. 3: Fortli-roo t of t he condensate propagatio n of t he re- 
duced densities ^ pi{t) (dark gray surface) and y/psii} (light 
gray surface). The non-adiabatic effects, seen as population 
escaping the main wave packet, become clear in this plot by 
presenting the fourth-root of the densities. The parameters 
and conditions (given in the text) are the same as for fig-EKa) 
and (b). 



W{x) is for the |1) atoms, and consequently the critical 
velocity is instead around 35 cm/s. Again, both these 
velocities could be increased by using stronger lasers. 
The forth-root of the reduced densities, i.e. -y/ pi (t) and 
psit), for the positively propagating case are displayed 
in fig. [3] (the figure only gives one out of the 500 quan- 
tum Monte-Carlo trajectories). By showing the fourth- 
root of the densities, the small non-adiabatic contribu- 
tions during the population transfers are seen. These ap- 
pear as a small amount of population being trapped in 
the first STIRAP interaction region. During the second 
STIRAP, information leaks the system due to coupling of 
the cavity to its environment. The resulting decoherence 
acts as a reduction of the system wave function. This 
decoherence-induced collapse actually increases the effi- 
ciency of the diode. From fig. [H it is especially clear 
that the non-adiabatic contributions established during 
the first STIRAP process are suppressed during the wave 
function reduction. Due to the non-adiabaticity, the mo- 
tional degrees-of-freedom of the condensate become entan- 
gled with the cavity field. However, once the field decays, 
the entanglement is lost and the motional wave packet of 
the condensate collapses at the same instant. This effect 
becomes more transparent at high velocities. For a con- 
densate with an initial velocity of 50 cm/s, 74 % of the 
population remains in the dark state after the first STT 
RAP, but during the decoherence process of the second 
STIRAP the system relaxes almost perfectly back to its 
target state |1). Thus, even for non-adiabatic evolution, 
the scheme works extremely well. 



It is known that nonlinearity can greatly suppress adi- 
abaticity in the STIRAP scheme [55], as well as in more 
general settings \27\ . For sufficiently strong non-linearity, 
e.g. increasing the atomic density or the scattering length 
Us, the adiabatic potentials/energies form a loop-structure 
which prevents full adiabaticity regardless of how long the 
characteristic time-scale is [5^. The breakdown of adi- 
abaticity originating from atom-atom interaction was re- 
cently demonstrated in a double- well BEC system [28] . 
We have numerically checked the diode efficiency by vary- 
ing the atom number TV (equivalent to keeping the density 
fixed while tuning the scattering length as). For the given 
value of Qs and within the range between 10 000 to 2 000 
000 atoms, we do not find substantial changes in the ef- 
ficiency. Thus, we conclude that within this range, and 
for the present parameters, there are no loop-structures 
formed. 

Conclusion. — We have generalized the atom diode 
of ref. [To] to apply for a condensate of ultacold atoms. 
The idea of [TU] relies on irreversible spontaneous emis- 
sion which in a condensate will cause heating as well as 
decoherence and thereby breakdown of condensation. The 
present scheme reckon upon coherent coupling of the con- 
densate with a cavity mode, which in return couples irre- 
versibly to a photon reservoir. Decoherence in the conden- 
sate is then induced via the cavity field decay and should 
consequently be considerably suppressed. More precisely, 
we demonstrated the method on a mean-field level of the 
BEC, which does not fully take undesired quantum noise 
of the true atomic many-body system into account. Nev- 
ertheless, we believe that in terms of atomic condensates 
the present method should be superior the one of ref. [10] , 
and photon reservoir noise has indeed been considered in 
the present model. 

Utilizing experimentally relevant parameters, it was 
particularly shown how almost perfect efficiency could be 
reached. The system also gave an interesting example of 
non-locality and how the reservoir-induced wave packet 
reduction improves the diode performance, e.g. non- 
adiabatic contributions are annihilated by the collapse. 
The model Hamiltonian is, of course, an idealized case, 
and we have not taken; atomic many-body effects, heat- 
ing due to thermal atoms, nor system fluctuations into 
account. On the other hand, the process time-scale is at 
the ms level, which is well below typical survival times 
of atomic condensates. It should be pointed out that the 
strength of the cavity decay k sets some conditions on 
repeatability of the diode; namely if a second condensate 
passes through the diode it is assumed that the cavity field 
has had enough time to relax to vacuum. With the k used 
in our examples this would certainly be true as the cavity 
field amplitude is approximately zero as soon as the con- 
densate has left the interaction region. Finally, the present 
scheme relies on the STIRAP process, but we note that 
possible other adiabatic processes could be employed, for 
example a Landau-Zcner transition j29| between the two 
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atomic states |1) and |3). However, such chirp-scheme 
woTiki typically require changing of the internal atomic 
energies, which is most likely more difficult to accomplish 
experimentally. Therefor, the present STIRAP proposal 
is probably more relevant for current experimental setups. 

* * * 
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